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ABSTRACT 

We analyse the QCD chiral phase transition in the nonhnear and hnear a- 
model. The strategy is the same in both cases. We fix the parameters of the 
effective meson theory at temperature T = and extrapolate the models to 
temperatures in the vicinity of the phase transition. The linear cj-model in 
SU{3) X SU{3) gives a crossover around Tc ~ 190 MeV. Around this temperature 
chiral SU{2) x SU{2) is almost restored. We also calculate meson masses as a 
function of temperature. 

1. Introduction 

Relativistic nucleus- nucleus collisions with cm-energies Ecm > 20 GeV/nucleon 
create large systems of sizes i? > 20 fm at freeze-out with > 10^ pions Jll.It is natural 
to try statistical methods to describe such hadronic fireballs. A good starting point may 
be to use equilibrium thermodynamics with pions when one is interested in the later 
stages of these collisions, where the temperature is below possible phase transitions. 
The low energy interaction of pions is fully determined by chiral symmetry [0,0. 
Below T K, 120 MeV this interaction can be parametrized by the nonlinear sigma model 
0(4) with one scalar (a) and three pseudoscalar vf-fields, which are constrained by the 
condition + vf^ = where f^ is the pion decay constant. Above this temperature 
region also heavier hadrons give a non-negligible contribution to the condensates and 
thermodynamic quantities. One way of including part of the heavier mesons is provided 
by the choice of SU{?>) x SU{?>) as chiral symmetry group rather than SU{2) x SU{2). 
The linear SU{?>) x SU{?>) sigma model includes a nonet of pseudoscalar (0~)-fields 
and a nonet of scalar (O^)-fields. The spontaneous breaking of the S'f/(3) x SU{'i) 
symmetry leads to massless (O^) Goldstone modes. Obviously a massless pseudoscalar 
octet does not provide an adequate approximation to the experimentally observed 
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meson spectrum. Therefore we include explicit symmetry breaking terms to account 
for the physical mass values of the octet-fields. A determinant term in the meson 
fields guarantees the correct mass splitting of the 77 — 77' masses which is due to the 
t/(l)-anomaly. It reflects the 't Hooft-determinant on the quark level. 

The experimental challenge is to measure the equation of state of pions from the 
inclusive pion spectra. The theoretical task is to calculate this equation of state. For 
this purpose we need reliable techniques to treat field theories at finite temperature. 
One can then extrapolate from the measured physics at T = to the yet unknown 
physics at high temperatures. A very accurate treatment of the soft modes with lowest 
mass is essential at low temperatures. We calculate the partition function Z in terms of 
a selfconsistent field which is chosen to extremize In Z. It gives effective masses to the 
meson fields. This saddle point approximation to the partition function corresponds 
to the leading order of a 1/iV expansion. 

The paper is organized as follows. In section 2 we discuss the nonlinear cr-model for 
SU{2) X SU{2). In section 3 we calculate the partition function in the linear cr-model 
for SU{3) X SU{3). Section 4 is devoted to a short discussion. 

2. The Nonlinear a-Model: SU{2) x SU{2) 

The partition function Z for the SU (2) x SU (2) nonlinear cr-model is given in terms 
of the 0(4)-multiplet (no, n) — (cr, yf) with = cr^ -|- 7?^ as: 

Z^l Vn{x) n Kn\x) - fl) exp{- dr d'x[^{d,nf - crio]} . (1) 

At zero temperature T := l/f3 = the parameters of the model are well known. The 
pion decay constant equals to 93 MeV. The classical vacuum expectation value of 
tiq is determined as < tzq >= /tt by minimizing the vacuum energy. Expanding the 
dependent field no = ^ — ft'^ to leading order in fp/ f^, one obtains the mass of the 
pion as ttt.^ = c/Ztt. 

The basic idea of our method is to eliminate the nonlinear constraint n"^ = a'^+n^ — 
fl by the introduction of an auxiliary field \[x) via an integral along the imaginary 
axis 

a+ioo ^ 

Z^ j VX{x) lvn{x)exp{- I dr J d^x[-{d^n)^ + X{n^ - f^) - crio]} . (2) 

a—ioo 

After shifting the zeroth component uq to uq 

~ c 
h = (no,n) = (rzo " ^'^) ' 

we obtain a Gaussian action for the 0(yN) multiplet field n, when we evaluate Eq.(2) 
in a saddle point approximation. The resulting partition function is given as 

Z^ j Vh{x) exp{- dr d^x[\{d,nf + XfP - Xfl - |^]} . (3) 
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Here we have dropped the A-integration and chosen X{x) = X — const. The optimal 
choice for A will be determined later, cf. eq. (11). 

Upon Gaussian integration over the four (A^ = 4) fi-fields we end up with a partition 
function of a free relativistic Bose gas with = 4 components and effective masses 

miff = 2A. (4) 

Z = eM-mUrimlff, A) + Uo{mlff, A) - Xf, - ^]} . (5) 
Here Ut denotes the contribution from thermal fluctuations 

Urimlff, A) = 4- ^ ln(l - e (6) 



and Uq the zero point energy 



We regularize the /^-integrations in Eqs. (7) and (8) with a cut-off A, since we do 
not believe our 7f-effective theory to be correct for momenta beyond A. At momenta 
A; > A the compositness of the pions manifests itself in resonance excitations and/or 
higher derivative couplings of the pion states, which are neglected. For the numerical 
calculations we take cut-off values A = 700 MeV, 800 MeV, 1000 McV. 

After regularization we adopt the following renormalization procedure. We define 
a renormalized potential at arbitrary T according to 



W^-{mlff,k) = UT{mlff,k) + U,{mlff,k) 



-[Uo{ml,K) + {mlff-mi) 



dmlff 



(8) 



The two subtractions guarantee the two renormalization conditions at T = 



m%f{X,) = ml (9) 
<nQ> = U . (10) 

It is well known that the nonlinear sigma model is not renormalizable in four di- 
mensions. Therefore higher order divergences can only be compensated by higher order 
derivative terms in the original action. The coefficients of these higher order terms have 
to be determined by experiment. We do not include such terms in contrast to ref. [2]. 
In section 3 we will extend the calculation to the hnear cr-model {SU{Z) x SU{Z)) 
which contains higher masses and strange mesons. 
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The thermodynamic observables at finite T are obtained from the partition function 
Z approximated as 



Z(A,T) = exp<^ -/?\/ 



f/-K^/A*),A)-A*/.'-^ 



where A*(r) extremizes \nZ aX a given temperature T 7^ 0. The saddle point equation 
for A* is solved numerically, since in the interesting temperature range the relevant 
parameter rrieff/T = v2A*/T can have values in the range < nieff/T < 00. 

Let us first study, how the order parameter of chiral symmetry breaking < no > 
behaves as a function of temperature. In Fig. 1 we present the result for = 

^^'^dc^^ I ^^"i^J""'' • ■'■^ quark language this ratio corresponds to the ratio of the quark 
condensate < qq{T) > at finite temperature over the quark condensate at T = 0, since 
the symmetry breaking term of the 0(4) Lagrangian Csb = cuq is identified with the 
symmetry breaking term £sb = —2mqq in the QCD-Lagrangian. We also show the 
result of the linear cr-model SU (3) x SU (3), which are presented in the next section, and 
the results of chiral perturbation theory 0. The result for< qq{T) > / < qq{T = 0) > 
is rather insensitive to the cut-off. It agrees well with the three loop calculation of 
ref. [2]. Chiral symmetry is only very gradually restored. At low temperature the 
TTTT-interaction is weak and < qq > does not change very much. 

3. The Linear cr-Model: SU{3) x SU{3) 

For a Euclidean metric the Lagrangian of the linear sigma-model is given as 

C = ^tra^$a^$+-^/i2tr$$+ + /i(tr$$+)' + /2tr ($$+)' 

+5((det<l> + det<l>+) -^oc^o-^sCTs, (12) 



where the (3 x 3) matrix field ^(x) is given in terms of Gell-Mann matrices A^ {i = 
0, . . . , 8) as 

1 s 

Here ai and 7i£ denote the nonets of scalar and pseudoscalar mesons, respectively. As 
order parameters for the chiral transition we choose the meson condensates < ctq > and 
< CTg >. The chiral symmetry of C is explicitly broken by the term (— eo^o — £8<^8), 
corresponding to the finite quark mass term 2mqqq + rrisSS on the quark level. The 
chiral limit is realized for vanishing external fields Eq and es- Note that the action 
S = J dPxdrC with C of Eq. (|I2[) may be regarded as an effective action for QCD, 
constructed in terms of an order parameter field $ for the chiral transition. It plays 
a similar role to Landau's free energy functional for a scalar order parameter field for 
investigating the phase structure of a $^-theoryH]. 
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The six unknown couplings of the sigma-model (Eq. ([T^)) {fi^, fi, f2, g,eo,e8) are 
assumed to be temperature independent and adjusted to the pseudoscalar masses at 
zero temperature. Further experimental input parameters are the pion decay constant 
/tt = 94 MeV and a high lying (0+) scalar mass rricr^ = 1.20 GeV (cf. Table 1). For 
the remaining scalar masses and the coupling constants we obtain the values of Table 
1. 



Table 1: Tree level parametrization of the SU{3) x SU{3) linear sigma model. 

Input 



[MeV] 


rriKlMeV] 


m^[MeV] 


rrin' [MeV] 


U [MeV] 


= ^h [MeV] 


130. 


491. 


544. 


1046. 


94. 


1200. 


Output 




/i 




9 [GeV] 


eo [GeV^] 


es [GeV^] 


0.059 


4.17 


4.48 


-1.8 


0.026 


-0.035 




niao [GeV] 


m^. [GeV] 


mf> [GeV] 








1.0116 


1.0312 


0.7495 









The interpretation of the observed scalar mesons is controversial. There are good 
reasons to interpret the (0"*") mesons at 980 MeV as meson bound states. The model 
underestimates the strange quark mass splitting in the scalar meson sector, the value 
for TTiK* comes out too small. 

The effective theory can be related to the underlying QCD Lagrangian by comparing 
the symmetry breaking terms in both Lagrangians and identifying terms with the 
same transformation behaviour under S'f/(3) x SU{?)). Taking expectation values in 
these equations we obtain the following relations between the light quark condensates, 
strange quark condensates and meson condensates 

£^o(o-o) = -^(2m + m,)(2(gg) + (ss)) 
2 

£8(fX8) = -- (m - m,) (2 (gg) - (ss)) . (14) 

Weusem= (m„ + mrf)/2 = (11.25 ±1.45) MeV and = (205 ±50) MeV for the light 
and strange quark masses at a scale A = 1 GeV P] . From the scalar meson condensates 
at T = 0, ((To) = 0.14 GeV and (ag) = -0.03 GeV we get 



{qq) = - (243.8 ± 60MeV)^ 
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{ss) = - (285.0 ± 30MeV)^ 



(15) 



in accordance with values from PCAC relations within the error bars. Since we 
treat the coefficients Eq, Es of < (Jq > and < as >, and m, of < gg > and < ss > as 
temperature independent, we will use Eqs. (|T^) for all temperatures to translate our 
results for meson condensates into quark condensates. 

We also check that the pseudoscalar meson mass squares, in particular and rnj^ 
are linear functions of the symmetry breaking parameters Eq,Es- Varying ^8 while 
keeping the other couplings fixed we can simulate the sigma model at unphysical meson 
masses. Since the current quark masses are assumed to depend linearly on Eq and Eg, 
an arbitrary meson mass set can be related to a mass point in the (m^^rf, ms)-plane by 
specifying the choice of {eq, Eg)- This may be useful in order to compare our results for 
meson (and quark) condensates with lattice simulations of the chiral transition. 

The thermodynamics of the linear sigma model is determined by the partition 
function with the Lagrangian of Eq. (p^) 

Z = Jv^expi^-J^^dTjd^xCmx,T))Y (16) 

We will treat Z again in a saddle point approximation. As mentioned above, the saddle 
point approximation amounts to the leading order of a 1/A^-expansion. In this model 
N = 2N] = 18. Note that C of Eq. (|12D would be 0(A^)-invariant, if = and 
g = 0. Our input parameters lead to non-vanishing values of /2 and g, therefore the 
0(A^)-symmetry is only approximately realized. 

We calculate the effective potential as a constrained free energy density Ueff{C,o, (.s), 
that is the free energy density of the system under the constraint that the average 
values of (Jq and ag take some prescribed values i^o and ^g- The values ^o^^^ and 
^8^;^ that minimize Ueff, give the physically relevant, temperature dependent vacuum 
expectation values, i.e. < cxo >= ^o™i„, < cts >= ^Smin- Hence we start with the 
background field ansatz 

O"0 = ^o + o-Q 

(TS = + (17) 

where a'^ and CTg denote the fluctuations around the classical background fields and 
^8- AH other field components are assumed to have zero vacuum expectation value, i.e. 
o"£ = 0"^ for ^ = 1, . . . , 7 and vr^ = vr^ for £ = 0, . . . , 8. The relation between the effective 
potential Ueff and Z is given by 

Z = J d^o J d^se-^^^^"'^^"'^'^ . (18) 

Next we insert the background field ansatz (0) in C and expand the Lagrangian 
in powers of $' = ^ J2£=oi'^£ + ^^^)'^^- The constant terms in $' lead to the classical 
part of the effective potential Udass- Linear terms in $^ vanish for all £ = 0, . . . , 8 due 
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to the 5-constraints in Eq. (|T8|). Quadratic terms in define the isospin multiplet 
masses mq, where Q = 1, . . . , 8 labels the multiplets. 

The cubic part in $' will be neglected, while the quartic term C^'^\^') is quadratized 
by introducing an auxiliary matrix field Y^{x, r). This is a matrix version of a Hubbard- 
Stratonovich transformation 0. 

In the saddle point approximation we drop / PS and use a SU (3)-symmetric ansatz 
with a diagonal matrix = diag(s, s, s). Hence the effect of the quadratization proce- 
dure is to induce an extra mass term (s + /Xq) and a contribution [/saddle to Ueff, which 
is independent of and ^g. This way we finally end up with the following expression 
for Z 



g 



where 



Xl = ml + fxl + s, (20) 

ip'q denotes Uq for Q = 1, ... ,4 and ttq for Q = 5, . . . , 8, g{Q) is the multiplicity 
of the isospin multiplet. We have g{l) = 3 for the pions, g{2) =4 for the kaons, 
g{3) = 1 = g{4:) for rj,-)]', respectively. Correspondingly, the multiplicities for the scalar 
nonets are g{5) = 3, g{6) = 4, g{7) = 1, g{8) = 1 for the oq, Kq, /q, /g-mesons. 

Thus we are left with an effectively free field theory. The only remnant of the 
interaction appears in the effective mass squared Xq via the auxiliary field s. 

The choice of a self-consistent effective meson mass squared has been pursued al- 
ready in Refs. ||^, Q . This is an essentially new ingredient compared to earlier calcula- 
tions of the chiral transition in the linear sigma model 0. The positive contribution 
of s to the effective mass extends the temperature region, where imaginary parts in 
the effective potential can be avoided. In general, imaginary parts are encountered, 
when the effective mass arguments of logarithmic terms become negative. They are 
an artifact of the perturbative evaluation of the effective potential and of no physical 
significance, as long as the volume is infinite. In our application the optimized choice 
for s will increase as function of temperature and lead to positive Xq over a wide range 
of parameters. 

Gaussian integration over the fiuctuating fields in Eq. (|l^) gives 

Z{^0, ^8, S) = exp {-PV [Ucl^s + [/saddle + 

4g*W)E/||>"(^^K + -J))l} (21) 

where 

u;'Q^e + X'q, (22) 
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and 

^ (2Tin/Tf (23) 

denote the Matsubara frequencies. In contrast to our former approach P] we keep all 
Matsubara frequencies and evaluate I]nez iii the standard way, see e.g. [|1^]. The result 
is 

Ziio.is-.s) = e-^^^-«-^«^^) (24) 

Uesiio,^?,'-, s) = f/class + f^addle + t^h (25) 

- ig««)/(05>n(l-.--) (26, 
Here we have indicated that Z and Ucs still depend explicitly on the auxiliary field 

s. 

The linear sigma model is a renormalizable theory, and in principle the zero point 
energy can be calculated and renormalized. We do, however, believe that this model 
is an effective description for QCD already at tree level. 

Now we are prepared to determine the temperature dependence of the order pa- 
rameters < ^0 > (T), < ^8 > (T) from the minima of Ues{^o,^8', s*). Thermodynamic 
quantities like energy densities, entropy densities and pressure can be derived from Z 
in the standard way, if Z is approximated as 

2 = g-/3V'(7eff(5o,C8;s*)_ ^27) 

For the parameters of Table |l] we vary the temperature and determine for each T 
the extremum of Ues as a function of ^8 and s. The extremum is a minimum with 
respect to and and a maximum with respect to s. For the search of the saddle point 
it is necessary to continue the effective potential into the region of complex effective 
masses Xq (cf. Ref. |Tl 



In Fig. 2 we show the variations of ^^^^^q and ^^^^'^q as a function of temperature 
obtained from < > (T) and < ,^8 > {T) with the help of Eq. (|np. We observe a 
gradual decrease of the light quark condensate, whereas the strange quark condensate 
stays almost constant. 

In our lowest order calculation the temperature dependence of these masses is de- 
termined by the temperature dependence of the condensates [cf. Fig. 3]. The masses 
mi and become degenerate after the crossover. The vr — K- splitting is increased 
rather than reduced. Accordingly the strange meson contribution to the energy density 
in this temperature region is reduced compared to the low- temperature hadron gas. 

In Fig. 4 we give the energy density u/T^, entropy density s/T^ and pressure p/T^- 
Sizeable contributions to u come mainly from 8 degrees of freedom, the pions, the kaons 
and the /q meson. 
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4. Discussion of the Results 

For low temperatures the physics of the nonlinear SU (2) x SU (2) and linear S't/(3) x 
SU (3) cr-model are identical. In Fig. 1 we show the light quark condensates calculated 
in both models. Above T ^ 120 MeV the extra degrees of freedom in the SU (3) x SU{3) 
calculation become important. At higher temperatures T ^ Tq ~ 190 MeV also 
the linear sigma model will certainly fail as an effective model for QCD due to the 
lack of quark-gluon degrees of freedom. Nevertheless it would be interesting to study, 
at what temperature the full SU{3) x SU{3) symmetry is restored. At very high 
temperatures the effective potential becomes proportional to J2q ^'^{Q)T'^, the linear 
terms proportional to ctq and as in the masses of O'^ and 0~ mesons cancel and 
temperature tries to fully restore the broken symmetry. 

Finally we remark that the crossover occurs around Tq ~ 190 MeV, which is rather 
close to the Hagedorn temperature Th{Th ~ 160 MeV). This may not be entirely 
accidental. In our model the 1/A^-expansion means a large number of flavours, since 
N = 2 ■ Nj. In order to keep QCD an asymptotically free theory also the number of 
colours Nc has to increase. Correspondingly our approximation is similar to Hagedorn's 
description of the hadron gas as a resonance gas. We expect that corrections from 
subleading terms in our l/iV/-expansion will implicitly amount to corrections also to 
the large Nc-limit. The chiral transition for unphysical values of strange quark and 
light quark masses will be investigated in a future publication |T^. 



Figure Captions 

Fig. 1: The order parameter for chiral symmetry breaking {qq{T)) / {qq{T = 0)) for 
two different cut-offs A = 700 MeV (solid line) and A = 1000 MeV (dashed 
line). The diamonds represent the result of the linear a- model in section 3. 
The crosses are the result of ref. 2. 

Fig. 2: Normalized light quark condensate (qq) / {qq)2^^Q vs temperature (full curve) 
and normalized strange quark condensate (ss) / {ss)rp^Q vs temperature (dash- 
ed curve). 

Fig. 3: Meson masses Xq, Q = 1, ... ,8 as a. function of temperature T. {X^^ = mag , 
X„^ = tuk* , A:^,,, = ruf,^ , X„^ = m/J. 

Fig. 4: Entropy density s/T^, normalized energy density m/T^ and pressure p/T^ vs 
temperature. Errors are only indicated for s/T^. 
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